Crossover and self-averaging in the two-dimensional site-diluted Ising model: 
Application of probability-changing cluster algorithm 
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Using the newly proposed probability-changing cluster (PCC) Monte Carlo algorithm, we sim- 
ulate the two-dimensional (2D) site-diluted Ising model. Since we can tune the critical point of 
each random sample automatically with the PCC algorithm, we succeed in studying the sample- 
dependent Tc{L) and the sample average of physical quantities at each Tc{L) systematically. Using 
the finite-size scaling (FSS) analysis for Tc{L), we discuss the importance of corrections to FSS both 
in the strong-dilution and weak-dilution regions. The critical phenomena of the 2D site-diluted Ising 
model are shown to be controlled by the pure fixed point. The crossover from the percolation fixed 
point to the pure Ising fixed point with the system size is explicitly demonstrated by the study of 
the Binder parameter. We also study the distribution of critical temperature Tc{L). Its variance 
shows the power-law L dependence, and the estimate of the exponent n is consistent with 

the prediction of Aharony and Harris [Phys. Rev. Lett. 77, 3700 (1996)]. Calculating the relative 
variance of critical magnetization at the sample-dependent Tc{L), we show that the 2D site-diluted 
Ising model exhibits weak self-averaging. 

PACS numbers: 75.50.Lk, 64.60.Ak, 05.10.Ln 



I. INTRODUCTION 



The critical behavior of random spin systems have been 
studied for more than decades |^ . The effect of random- 
ness has been attracting much attention because real ma- 
terials contain impurities. In a pioneering work, Harris 
studied the problem whether the randomness changes 
the critical behavior of the pure system in the case of 
quenched random systems. The so-called Harris crite- 
rion 1^ states that even a weak randomness changes the 
critical behavior if the specific heat exponent apuic of the 
pure system is positive; the random system belongs to the 
different universality class from that of the pure system. 
Then, the randomness becomes relevant in the renormal- 
ization group (RG) terminology; the critical behavior is 
controlled by the random fixed point. Earlier works on 
random Ising models |^,^ deserve to be mentioned. 

The two-dimensional (2D) Ising model is a marginal 
case. The specific heat diverges logarithmically at the 
critical point, that is, apurc = 0. We cannot tell whether 
the randomness is relevant or not from the Harris crite- 
rion. Using the quantum-field theory, Dotsenko and Dot- 
senko |^ calculated the correlation length of the 2D di- 
luted Ising system. They showed that the randomness is 
irrelevant but there appear logarithmic corrections. The 
expression for the correlation length in the weak-dilution 
region is given by 



(l + Aln(l/|^|))V^ 
\t\ 



t={T- T,)/T„ 



(1) 



(l + Aln(l/|t|))i/i6' 



(2) 



where A(<C 1) is the strength of randomness. Shalaev, 
Shankar, and Ludwig proceeded with the calcula- 

tion of magnetization, 



which also shows the logarithmic corrections to the pure 
Ising model without any change in critical exponents. 

The crossover is an interesting subject especially for 
random spin systems. The critical behavior is controlled 
by a relevant fixed point. However, it is influenced by 
irrelevant fixed points. Thus, as the system size becomes 
large, we observe the crossover behavior at the critical 
region between the relevant fixed point and irrelevant 
fixed points. In the case of diluted spin models, the mag- 
netic order disappears at the percolation threshold even 
at r = 0. Then, the crossover between the pure fixed 
point, the random fixed point (if there exists), and the 
percolation fixed point is the subject of concern ||]. 

The problem of self-averaging is also of current inter- 
est for random spin systems |10 1^ because each sample 
has a different random configuration. The system is said 
to exhibit self-averaging if the relative variance of the 
thermal average of a quantity goes to zero as the system 
size becomes infinite. Using the RG, Aharony and Har- 
ris (AH) discussed the L dependence of the relative 
cumulants of singular quantities, where L is the linear 
system size; they discussed the self-averaging property 
of the random system, which depends on whether the 
randomness is relevant or not. 

The Monte Garlo simulation is a powerful tool to study 
difficult problems such as random spin systems. However, 
the simulation method sometimes suffers from the prob- 
lems of slow dynamics. Several new algorithms for the 
Monte Carlo simulation have been proposed to overcome 
such difficulties. Gluster algorithms (l2|l2] are examples 
of such efforts. The histogram method |14] enables us to 
calculate physical quantities for different parameters with 
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a single simulation by using the reweighting technique. 

Quite recently the present authors have proposed 
a new effective cluster algorithm, which is called the 
probability-changing cluster (PCC) algorithm jl^, of 
tuning the critical point automatically. The invaded 
cluster algorithm [|l6| was also proposed to determine 
the critical point automatically, but its ensemble is not 
necessarily clear. In contrast, with the PCC algorithm 
we approach the canonical ensemble asymptotically; we 
can use the finite-size scaling (FSS) analysis for physical 
quantities near the critical point. The PCC algorithm 
is quite useful for studying the random spin systems, 
where the distribution of the critical temperature, Tc, 
due to the randomness is important, because we can tune 
the critical point of each random sample automatically. 
Wiseman and Domany [ pH pointed out the importance 
of determining the sample-dependent Tc{L) in the simu- 
lational study of random spin systems. They applied a 
reweighting technique to search critical points of random 
samples, but the iteration process was needed and there 
was a difficulty in tuning Tc{L) for some exceptional sam- 
ples. The importance of the random average a/ter finding 
the critical point of each sample was also pointed out by 
Bernardet et al. p7| ]. 

In this paper, we study the 2D site-diluted Ising model 
using the PCC algorithm. We focus on the crossover phe- 
nomena and the sample dependence of physical quanti- 
ties. The rest of the paper is organized as follows. In Sec. 
II, we describe the model and the simulation method, the 
PCC algorithm. In Sec. Ill, we address the FSS analysis 
of the data and the crossover between fixed points, pay- 
ing attention to the corrections to FSS. In Sec. IV, we 
study the variance of the quantities, and discuss the self- 
averaging property of the 2D site-diluted Ising model. 
The summary and discussion are given in Sec. V. 



of the Swendsen-Wang cluster algorithm Jl^ , where the 
Kasteleyn-Fortuin (KF) representation |26|] of the Ising 
(Potts) model is used. To form a cluster, parallel spins 
are connected with the probability p^^ = 1 — g^-'^J/f^BT _ 
For the diluted model, of course, only pairs of spin sites 
are connected. In the PCC algorithm we change the 
probability depending on the observation whether 
clusters are percolating or not percolating. A simple neg- 
ative feedback mechanism together with the FSS prop- 
erty of the percolation leads to the determination of the 
critical point. As Ap^^, the amount of the change of p^^, 
becomes small, the distribution of becomes a sharp 
Gaussian distribution around the mean value p^^(i), 
which results in the determination of the critical tem- 
perature Tc{L). We approach the canonical ensemble in 
this limit. For the more detailed description of the PCC 
algorithm, see Ref. |l5|] . 

We have made simulations for the 2D site-diluted Ising 
model on the square lattice with the system sizes L = 16, 
32, 64, 128 and 256. For the spin concentration p, we have 
picked up p = 0.7 and p = 0.9. We deal with the grand- 
canonical ensemble of samples in which the occupation 
of each site is determined independently with probabil- 
ity p. The percolation threshold for the site percolation 
of the square lattice is known to be 0.592746 p^. The 
final value of Ap^^ has been chosen as Ap^^ = a/L, 
where a — 6.25 x 10~^ for the system size L. There are 
several choices for the criterion to determine percolating. 
We have employed the topological rule in the present 
study. The topological rule is that some cluster winds 
around the system in at least one of the D directions 
in _D-dimensional systems. The random average is taken 
over 2,000 samples for L = 16 to 128, and 1,000 samples 
for L = 256. For each random sample we have made 
10,000 Monte Carlo sweeps to take the thermal average. 



II. MODEL AND SIMULATION METHOD 

We are concerned with the site-diluted Ising model 
whose Hamiltonian is given by 

TL ^ -j'^tiejSiSj. (3) 

Here Si is the Ising spin on the lattice site i, and is 
a random variable that takes 1 (spin) or (vacancy). 
The summation is taken over the nearest-neighbor pairs 
(i, j). The concentration of the spin will be denoted by p. 
The 2D site-diluted Ising model was already studied by 
the Monte Carlo methods |Q , and the 2D random-bond 
Ising model was investigated by the Monte Carlo sim- 
ulations ]l7| , p^| -p3[ , the transfer-matrix calculation [ p4[ , 
and the high-temperature expansion [25|] . In the present 
study we give special attention to the distribution of 
physical quantities using the PCC algorithm. 

Here we briefly describe the idea of the PCC algo- 
rithm |l^. The PCC algorithm is an extended version 



III. FINITE-SIZE SCALING ANALYSIS AND 
CROSSOVER 

Let us start with the size dependence of the critical 
temperature. We plot \Tc{L)] as a function of 1/L for 
both p = 0.7 and p = 0.9 in Fig. |[ where the brackets 
[• • ■] represent the random sample average. From now on, 
we represent the temperature in units of J /ks- The error 
bars are within the size of marks. 

We employ the FSS analysis for [rc(i)]. According to 
the theory of FSS ||2^, if a quantity Q has a singularity 
of the form Q{t) ^ near the criticality t = 0, the corre- 
sponding quantity Q{L, t) for the finite system with the 
linear size L has a scaling form 

Q{L,t)^L-^'^f{tL'/''), (4) 

where v is the correlation- length exponent. Then, if the 
corrections are negligible, the critical temperature for 
the infinite system can be estimated through the relation 

[T,{L)]=T,^AL-^'r (5) 
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But the corrections to FSS due to irrelevant fixed points 
are important especially in the strong-dilution region. 
Then, we use the relation 

[Te(L)] - + + BiL-'^) (6) 

to determine Tc and v for p = 0.7. Using the non-linear 
fitting, we estimate Tc as 1.0712(5) and \/v as 0.96(1) 
for p = 0.7. Here, the number in the parentheses de- 
notes the uncertainty in the last digits. The correction- 
to-scaling exponent is estimated as 0.63(20). In the 
weak-dilution region [p = 0.9), the corrections to FSS 
given in Eq. ^ are rather small. However, there ap- 
pear small logarithmic corrections, as was pointed out by 
Dotsenko and Dotsenko around the pure Ising fixed 
point; for finite systems the corrections to FSS due to 
the logarithmic term in Eq. (|l|) may be ascribed to the 
corrections to the critical exponent. Thus, we use the 
equation 

[Tc{L)] =Tc + aiL-i/'^+^'iC'^^)"'' (7) 

for the analysis of the data for p = 0.9. A similar analysis 
was employed by Talapov et al. [pT| with uj = 1. Here, 
we regard w as a free parameter to include higher-order 
contributions. Our estimate of Tc with the non-linear fit- 
ting is 1.9022(6) and that of l/u is 1.00(1) for p = 0.9. 
The solid curves in Fig. |l| are the best-fitted curves using 
Eq. (H) and Eq. (0) for p = 0.7 and p = 0.9, respectively. 
Since the critical exponent 1/v for the 2D Ising model 
is 1, our estimates of 1/v for the 2D site-diluted Ising 
model show that the critical phenomena are controlled 
by the pure Ising fixed point. 
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FIG. 1. Size dependence of critical temperatures \Tc{L)] 
(in units of J Iks) for p = 0.7 and 0.9. The system sizes L are 
16, 32, 64, 128, and 256. The solid curves are the best-fitted 
curves using Eqs. (^) and (Q) for p = 0.7 and 0.9, respectively. 

Next turn to the average value of magnetization at 
the critical temperature, [(|mc(L)|)]. Here, the angular 
brackets (• • •) represent the thermal average. We have 
measured the magnetization at the sample-dependent 
critical temperature Tc(L) for each sample; then we have 



taken the sample average. The efficiency of this proce- 
dure for high-precision analysis was pointed out by Wise- 
man and Domany ||l^. The log-log plot of [(|mc(i)|)] 
versus L for p — 0.7 and 0.9 is given in Fig. ^ The error 
bars are again within the size of marks. 
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FIG. 2. Logarithmic plot of [(j"ic(I/)|)] versus L for p = 0.7 
and 0.9. The system sizes L are 16, 32, 64, 128, and 256. 
The solid curves are computed by the non-linear fitting using 
Eqs. (^ and for p = 0.7 and 0.9, respectively. 

The FSS relation including the corrections to scaling 
for the magnetization may be described as 

[(|m,(L)|)]-A2L-'3/'^(l + B2L-") (8) 

for the strong-dilution region. We can estimate the expo- 
nent P/vhy using Eq. (||). We should note that we do not 
need the information of Tc for the estimate of (3 /v. We 
estimate the critical exponent fS/f as 0.123(2) forp = 0.7 
using the non- linear fitting. The estimated value of (i/v is 
consistent with that of the pure Ising model, 1/8=0.125. 
The estimate of the corrections-to-scaling exponent is 
0.66(4), which is the same as that for the case of [Tc(L)]. 
In the weak-dilution region, we take account of the loga- 
rithmic corrections For finite systems, we use the 
fitting function 

[(|m,(L)|)] =a2i-^/'^+''^('"^)"", (9) 

which is similar to Eq. (|^). The non- linear fitting based 
on Eq. (^ have led to the estimate of /S/iy as 0.125(1) for 
p = 0.9. 

For both cases of [Tc{L)] and [{\mc{L)\)], the FSS 
analysis taking account of the corrections to scaling has 
yielded the critical exponents of the pure Ising model, 
which shows that the critical phenomena of the 2D site- 
diluted Ising model are controlled by the pure fixed point. 
In the strong-dilution region, the corrections due to the 
irrelevant fixed point are important, whereas in the weak- 
dilution region, the logarithmic corrections around the 
pure Ising fixed point are dominant. 

The effects of irrelevant fixed points also appear as the 
crossover phenomena. In the present case, the percola- 
tion fixed point is such a fixed point. In order to see the 
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crossover explicitly, we study the size dependence of the 
Binder parameter ||29| defined as 



{m?) 



(10) 



which is often used in the analysis of the FSS. Since the 
prefactors of the L dependence in Eq. are canceled 
out, one may determine the critical point from the cross- 
ing point of the data of temperature dependence for dif- 
ferent sizes as far as the correction to FSS are negligible. 
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FIG. 3. Sample average of the Binder parameter [pj for 
p = 0.7, 0.9 and 1. The Binder parameter for the percolation 
problem is also shown. Solid curves are guiding the eye. 

Here, we investigate the sample average of the Binder 
parameter at each Tc{L), [gd- We plot [gd for p — 0.7, 
0.9 and 1 in Fig. ^, where p = 1 is nothing but the pure 
Ising model. The error bars are within the size of marks. 
As for the size dependence, the log L contributions may 
appear, although the explicit form of L dependence is not 
clear; we plot [gd as a function of l/log^QL in Fig. ||, 
and the solid curves are guiding the eye. We also show 
[gd for the percolation problem at the sample-dependent 
percolation threshold, Pc{L), in Fig. ^. We define the 
magnetization for the geometric percolation problem by 
assigning the Ising spins on each cluster randomly with 
the probability of 1/2. The values of [gd extrapolated 
to L = oo for p = 1 and percolation are shown by the 
arrows. We see from Fig. |3|that the Binder parameter of 
the dilute Ising model approaches the value of the pure 
Ising model as the system size becomes large. It is the 
manifestation of the crossover from the percolation fixed 
point to the pure Ising fixed point with the system size. 

The fraction of lattice sites in the largest cluster, c, 
plays a role of order parameter in the percolation prob- 
lem. The moment ratio of c, (c)^ / (c^), has the same FSS 
property as the Binder parameter because the prefactors 
of the L dependence are canceled out |^ . 

We plot the sample average [(c)^ / (c^)] at the criti- 
cal temperature of each sample for p = 0.7, 0.9 and 1 
in Fig. |[ For each sample, we have measured c at each 



Tc(L), and have taken average over random samples. We 
also give the data for the percolation problem. The error 
bars are within the size of marks. We show the values 
extrapolated to L = oo for p — I and percolation by the 
arrows. We find that the value of [(c)^ / (c^)] of the dilute 
Ising model approaches the value of the pure Ising model 
as the system size becomes large, which is the same as 
the case of the Binder parameter, [gd- 
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FIG. 4. Sample average of the moment ratio for the fraction 
of lattice sites in the largest cluster [(c)^ / (c^)] for p = 0.7, 
0.9 and 1. The data for the percolation problem are also 
shown. Solid curves are guiding the eye. 



IV. SELF-AVERAGING 

In this section, we study the distribution of physical 
quantities. We first deal with the variance of the critical 
temperature 



(Ar,(L))2 = [T,(L)2] - [T,(L)]^ 



(11) 



Using the RG argument, AH discussed the power- 
law size dependence, that is, (ATc(L))^ oc L~" for large 
L. According to AH, the exponent n depends on whether 
the random system is controlled by the pure fixed point 
or the random fixed point, that is, 



D (for pure fixed point) 

2/i'random (for random fixed point) 



(12) 



where D is the spatial dimension and ^'random is the 
correlation-length exponent for the random fixed point. 

We plot (Arc(i))^ for p = 0.7 and 0.9 in Fig. |. The 
variance {ATc{L))'^ becomes smaller as the system size 
becomes larger. Since the logarithmic scales are used in 
Fig. ^, we find the power-law size dependence from the 
linearity of the data; we can estimate the exponent n from 
the slopes of lines. Our estimates of n using least squares 
method are n = 1.95(3) for p = 0.7 and n — 2.02(3) for 
p — 0.9, respectively. These values are consistent with 
the prediction of AH [|ol, that is, n = D = 2. 
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FIG. 5. Size dependence of variance of critical temperature 
{ATc{L)f, given by Eq. Q. The slopes of solid lines give 
the exponent n for L""; n = 1.95(3) and 2.02(3) for p = 0.7 
and 0.9, respectively. 
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FIG. 6. Relative variance of critical magnetization at the 
sample-dependent Tc{L), R,n. Solid curves are guiding the 
eye. 



We treat the relative variance for the thermal average 
of physical quantity {X), 



Rx{L) 



[{Xf] [{X)] 

[{xw 



(13) 



when we discuss self-averaging. The system is said to 
exhibit self-averaging if Rx{L) ^ as L — ^> oo. AH ||l^] 
predicted that the size dependence of Rx {L) for the ran- 
dom system depends on whether the system is controlled 
by the random fixed point or the pure fixed point; that 



IS, 



Rx{L) cx 



(for pure fixed point) 



const, (t^ 0) (for random fixed point) 



(14) 



for L ^ oo. In the case of the random fixed point, 
the random system has no self-averaging. On the other 
hand, the system exhibits weak self-averaging in the case 
of the pure fixed point. Since the 2D Ising model is a 
marginal case, it is interesting to study the L dependence 
of Rx{L). 

We plot Rx (L) for the magnetization m at the sample- 
dependent Tc{L), Rrm as a function of L in Fig. ^. The 
absolute value of Rm is small for p = 0.9. It has larger 
value for smaller L for p = 0.7, but it becomes smaller 
for larger L. Since there may appear the logL contribu- 
tions, we have chosen 1/ logj^Q L for the horizontal axis as 
in Figs. ^ and ^, but the explicit form of size dependence 
is not clear. Anyway, from Fig. ^ we realize Rm — > 
as L oo; the 2D diluted Ising model exhibits weak 
self-averaging. In other words, the thermal average of 
magnetization does not depend on samples as the sys- 
tem size becomes infinite. 



V. SUMMARY AND DISCUSSIONS 

To summarize, we have investigated the 2D site-diluted 
Ising model using the PCC algorithm. Since we can tune 
the critical point of each random sample automatically, 
we have succeeded in studying the sample-dependent 
Tc{L) and the sample average of physical quantities at 
each Tc{L) systematically. 

We have used the FSS analysis for Tc{L). We have 
shown the importance of corrections to FSS due to irrel- 
evant fixed points in the strong-dilution region {p = 0.7); 
on the other hand, the logarithmic corrections are rele- 
vant in the weak-dilution region (p = 0.9). The critical 
phenomena of the 2D site-diluted Ising model are shown 
to be controlled by the pure fixed point with the loga- 
rithmic corrections. The crossover from the percolation 
fixed point to the pure Ising fixed point with the system 
size has been explicitly demonstrated by the study of the 
Binder parameter; it reflects the flow of renormalization. 

For the weak-dilution region, we have shown that the 
logarithmic corrections are important, which are compat- 
ible with the previous studies |l^ . We have also 
shown the crossover to the pure fixed point. It is opposed 
to the weak-universality scenario, that is, the exponent 
V is concentration-dependent |22 31]. No crossover was 
found in a recent work by Kim 23]. Here, we make a 
comment on the model. In the present study, we have 
studied the diluted Ising model. Ballesteros et al. jlj] 
studied the same model. The ground states are partial 
ferromagnetic states for dilution, and the magnetic or- 
der continuously disappears at the percolation threshold. 
The random bond mixture model of Ji and J2 couplings 
with 50 % concentration is also used ] p7pS| -j2^ because 
the critical temperatures are exactly known by the du- 
ality relation. The ground states of the latter model 
are always the complete ferromagnetic state except for 
J1I J\ = 0. In the weak randomness region, the correc- 
tions due to the randomness of the two models seem to 
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be the same; but the crossover behavior from the perco- 
lation point may be different. The careful study of the 
crossover for both models is still desired. Moreover, the 
replica symmetry breaking effects on the critical behavior 
of weakly disordered systems were also discussed 32 -3^] 
in relation to the Griffiths singularity The check of 
this problem will be left to a future study. 

We have studied the distribution of critical tempera- 
ture Tc{L). Its variance shows the power-law L depen- 
dence, and the estimate of the exponent, n ~ 2, is 
consistent with the prediction of AH, n ^ D. We have 
also calculated the relative variance of critical magneti- 
zation at the sample-dependent Tc{L), Rm- It becomes 
asymptotically close to zero as L becomes larger even for 
the case of p = 0.7. Thus, the 2D site-diluted Ising model 
is controlled by the pure Ising fixed point and exhibits 
weak self-averaging. 

In this paper we have studied the 2D diluted Ising 
model, where the pure Ising (relevant) fixed point and 
the percolation (irrelevant) fixed point are considered. 
For the three-dimensional (3D) diluted Ising model 
0,111- 37 1, the randomness becomes relevant because 



Opurc > 0; then there are three fixed points to be con- 
sidered, that is, the random fixed point, the pure Ising 
fixed point, and the percolation fixed point. The system- 
atic study of the 3D diluted Ising model using the PCC 
algorithm is quite interesting, which is now in progress. 
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